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Abstract 

An  algorithm  for  surface  domain  decomposition 
and  grid  generation  for  overset  grids  is  described.  A 
complex  surface  domain  is  covered  by  two  types  of 
grids:  seam  grids  and  block  grids.  Seam  grids  are  grids 
used  to  wrap  around  surface  crease  lines  and  regions 
of  high  surface  curvature.  Block  grids  are  grids  that 
cover  the  remaining  surface  regions  not  occupied  by 
the  seam  grids.  The  seam  and  block  grids  form  a  set  of 
overlapping  surface  grids  that  cover  the  entire  surface 
geometry.  Given  a  set  of  seam  grids,  an  automated 
method  for  generating  the  block  grids  is  presented  in 
this  paper.  Examples  of  seam  and  block  grids  are  given 
for  several  configurations  including  the  V-22  tiltrotor 
fuselage  and  the  X-CRV  Crew  Return  Vehicle. 

1.  Introduction 

The  generation  of  high  quality  grids  is  typically 
the  most  time  consuming  step  in  the  numerical  compu¬ 
tation  of  aerodynamic  flow  around  complex  configura¬ 
tions  using  structured  grids.  Amongst  the  structured 
grid  methods,  the  overset  grid  approach1  has  been  suc¬ 
cessfully  used  in  a  variety  of  applications.2-8  However, 
existing  grid  generation  software  requires  significant 
effort  from  the  user  to  generate  high  fidelity  surface 
grids  for  such  complex  cases.  This  is  primarily  due  to 
two  reasons.  The  first  reason  is  that  there  is  no  ro¬ 
bust  algorithm  or  software  available  to  decompose  the 
surface  domain  into  appropriate  sub-domains  for  sur¬ 
face  grid  generation.  Hence  this  step  depends  heavily 
on  the  experience  of  the  user  and  requires  intensive 
manual  effort.  The  second  reason  is  that  there  is  vir¬ 
tually  no  software  available  that  is  targeted  for  overset 
surface  grid  generation  where  a  marching  scheme9  is 
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frequently  more  convenient  than  an  algebraic  interpo¬ 
lation  or  elliptic  approach.  Since  the  marching  scheme 
described  in  Ref.  [9]  is  not  available  in  any  grid  gener¬ 
ation  system  with  CAD  capabilities,  overset  grid  users 
have  been  forced  to  use  software  designed  for  patched 
grids  such  as  GRIDGEN,10  ICEMCFD,11  and  NGP.12 

In  this  paper,  the  focus  will  be  on  methods  to 
improve  the  two  problem  areas  described  above:  sur¬ 
face  domain  decomposition  and  surface  grid  genera¬ 
tion.  There  have  been  several  efforts  to  automate  do¬ 
main  decomposition  for  the  patched  grid  approach.13,14 
These  require  the  decomposition  of  the  volume  space 
around  the  object  as  well  as  the  surface  of  the  object. 
This  strategy  is  not  ideally  suited  for  an  overlapped 
grid  approach.  Domain  decomposition  requirements 
for  overset  grids  can  be  reduced  to  a  surface  geometry 
decomposition  problem  where  the  surface  grids  are  al¬ 
lowed  to  overlap.  Once  the  surface  domain  is  decom¬ 
posed,  and  the  overlapping  surface  grids  generated,9 
body  conforming  volume  grids  can  be  produced  hyper- 
bolically  from  the  surface  grids  with  little  effort.15,16 
These  volume  grids  will  naturally  overlap  if  the  march¬ 
ing  distance  is  not  too  large.  The  rest  of  the  3-D  space 
can  then  be  filled  by  layers  of  Cartesian  box  grids.17’18 

Given  a  description  of  the  surface  geometry  (ref¬ 
erence  surface),  the  algorithm  for  surface  domain  de¬ 
composition  and  surface  grid  generation  presented  in 
this  paper  consists  of  two  steps:  seam  surface  grid  gen¬ 
eration  and  block  surface  grid  generation.  Seam  grid 
generation  is  primarily  a  manual  process  currently  and 
is  described  briefly  in  section  2.  Methods  for  perform¬ 
ing  block  grid  generation  automatically  are  described 
in  detail  in  section  3.  The  V-22  tiltrotor  fuselage  will 
be  used  as  an  illustrative  example  in  these  sections.  A 
panel  network  definition  of  the  V-22  fuselage  is  shown 
in  Figure  1.  Further  examples  of  seam  and  block  sur¬ 
face  grids  are  given  in  section  4. 
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Figure  1.  Surface  definition  of  V-22  tiltrotor  fuselage 
in  the  form  of  multiple  panel  networks. 

2.  Seam  Surface  Grid  Generation 

Seam  surface  grids  are  grids  bounded  by  one  or 
more  special  control  curves  (seam  curves)  of  the  sur¬ 
face  geometry.  The  generation  of  seam  surface  grids 
can  be  divided  into  four  steps.  Since  seam  grid  gen¬ 
eration  is  not  the  focus  of  this  paper,  each  step  is  de¬ 
scribed  only  briefly  below. 

(1)  Identification  of  seam  curves. 

The  seam  curves  are  special  control  curves  of  the 
surface  geometry.  These  control  curves  can  be  surface 
crease;  lines,  curves  along  high  curvature  contours  of 
the  surface,  intersection  curves  between  different  com¬ 
ponents,  and  other  special  curves.  Currently,  this  step 
is  performed  manually  but  there  is  good  potential  for 
automation  in  the  future. 

(2)  Identification  of  seam  surface  grid  domains. 

After  the  seam  curves  are  identified,  the  surface 
domains  that  are  to  be  bounded  by  one  or  more  of 
the  seam  curves  have  to  be  determined.  This  may  in¬ 
volve  splitting  or  concatenating  parts  of  different  seam 
curves.  Currently  this  step  is  performed  by  the  user;  it 
depends  heavily  on  the  user’s  experience,  and  is  per¬ 
haps  the  most  difficult  step  to  automate.  However, 
the  availability  of  a  graphical  user  interface  (GUI)  can 
significantly  reduce  the  human  effort  required. 


(3)  Distribution  of  grid  points  on  the  seam  curves. 

The  distribution  of  grid  points  on  the  seam  curves 
can  be  specified  directly  by  the  user  or  determined  au¬ 
tomatically  based  on  requirements  of  global  resolution 
and  local  geometric  constraints.  The  user  may  impose 
a  global  maximum  grid  spacing,  while  local  geometric 
features  may  demand  a  smaller  grid  spacing. 

(4)  Generation  of  seam  surface  grids. 

Finally,  the  seam  surface  grids  are  generated  from 
one  or  more  bounding  seam  curves.  For  a  given  seam 
surface  grid,  hyperbolic  or  algebraic  marching  schemes 
can  be  used  if  there  is  only  one  bounding  seam  curve. 
If  multiple  bounding  seam  curves  are  needed,  algebraic 
interpolation  schemes  or  elliptic  methods  can  be  used. 


Figure  2.  Seam  curves  and  seam  surface  grids  for  the 
V-22  tiltrotor  fuselage.  The  seam  curves  are  marked 
by  thick  lines. 

The  seam  curves  and  seam  surface  grids  for  the 
V-22  fuselage  example  are  shown  in  Figure  2.  A  seam 
grid  is  introduced  at  the  nose  to  avoid  the  need  for 
a  singular  axis  topology  for  a  surface  grid  at  this  lo¬ 
cation.  Other  seam  grids  are  introduced  at  the  crease 
line  on  the  fuselage,  the  wing/fuselage  intersection  line 
and  the  fuselage  trailing  edge.  All  grids  are  generated 
hyperbolically9  by  marching  from  the  seam  curves. 

3.  Block  Surface  Grid  Generation 

The  algorithm  used  to  automatically  construct 
block  grids  is  described  in  this  section.  It  is  assumed 
that  the  surface  geometry  definition  (reference  sur¬ 
face)  and  a  set  of  seam  surface  grids  are  supplied. 
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The  reference  surface  consists  of  one  or  more  refer¬ 
ence  patches  where  each  patch  can  be  a  Non-Uniform 
Rational  B-Spline  surface  (NURBS  surface),  a  panel 
network  (rectangular  array  of  surface  points),  or  some 
other  geometric  entity.  The  principles  presented  below 
are  independent  of  the  form  of  these  patches. 

A  software  program  called  SBLOCK  (stands  for 
Surface  BLOCKing)  has  been  developed  that  uses  mul¬ 
tiple  panel  networks  as  the  reference  surface  definition. 
The  panel  networks  are  assumed  to  be  abutting  but 
small  gaps  are  tolerated  between  different  networks. 
Point  to  point  match  between  adjacent  panel  networks 
is  not  required.  The  generation  of  block  surface  grids 
can  be  divided  into  six  steps.  Each  step  is  described  in 
more  detail  in  the  following  subsections.  Treatments 
for  other  types  of  reference  surface  definitions  such  as 
NURBS  surfaces  can  follow  similar  steps.  In  these 
cases,  the  reference  surface  is  first  converted  to  a  mul¬ 
tiple  panel  network  representation.  Then  the  steps  be¬ 
low  are  performed  followed  by  one  additional  final  step 
where  the  block  surface  grids  produced  are  projected 
back  onto  the  original  reference  surface  definition. 

3.1.  Creation  of  Marker  Cells 

The  first  step  in  block  surface  grid  generation  is 
to  create  a  set  of  marker  cells  on  the  reference  sur¬ 
face.  Each  marker  cell  is  just  a  small  region  of  the 
surface  domain  created  by  fracturing  the  surface  do¬ 
main  in  some  convenient  way.  Groups  of  these  cells 
can  be  flagged  to  mark  regions  where  block  surface 
grids  should  be  generated  as  described  in  section  3.2. 

Whether  the  reference  surface  is  prescribed  by 
patches  of  NURBS  surfaces,  panel  networks  or  some 
other  entities,  it  is  assumed  that  a  rectangular  ar¬ 
ray  of  points  can  be  created  to  represent  each  patch. 
For  patches  that  are  not  four-sided,  these  can  be  sub¬ 
divided  into  four-sided  sub-patches  where  a  rectangu¬ 
lar  array  of  points  can  be  placed  on  each  sub-patch. 
These  points  form  the  corners  of  a  rectangular  array  of 
quadrilateral  marker  cells  on  each  patch.  The  marker 
cells  can  be  made  to  be  uniformly  spaced  in  physi¬ 
cal  space  or  follow  some  other  constraints.  In  the 
SBLOCK  code,  points  on  each  panel  network  of  the 
reference  surface  are  redistributed  to  create  a  set  of 
uniform  marker  cells  such  that  the  size  of  each  marker 
cell  is  no  larger  than  a  user-prescribed  arc  length  A Sg. 


3.2.  Identification  of  Block  Cells 

The  next  step  is  to  identify  a  set  of  block  cells 
from  which  block  surface  grids  are  to  be  generated. 
Any  marker  cell  that  is  not  completely  enclosed  by  a 
seam  surface  grid  is  flagged  as  a  block  cell.  The  iden¬ 
tification  of  block  cells  is  accomplished  by  projecting 
the  vertices  of  all  the  marker  cells  onto  the  set  of  seam 
surface  grids.  Any  vertex  should  lie  either  inside  (in¬ 
cluding  the  boundary)  or  outside  of  a  seam  surface 
grid.  A  marker  cell  with  one  or  more  vertices  lying 
outside  of  all  the  seam  surface  grids  is  identified  as  a 
block  cell.  This  procedure  guarantees  overlap  between 
the  block  cells  and  the  seam  surface  grids.  The  block 
cells  identified  for  the  V-22  fuselage  example  are  shown 
in  Figure  3. 


Figure  3.  Block  cells  identified  for  the  V-22  tiltrotor 
fuselage. 

The  critical  element  in  this  step  (and  in  a  num¬ 
ber  of  steps  discussed  in  subsequent  sub-sections)  is  a 
robust  and  accurate  projection  routine.  Given  a  vertex 
of  a  marker  cell  (target  point),  the  set  of  seam  surface 
grids  that  may  contain  the  target  point  is  first  identi¬ 
fied  via  the  bounding  boxes  of  the  seam  surface  grids. 
Next,  the  seam  surface  grid  that  contains  the  closest 
point  to  the  target  point  is  selected.  A  Newton  stencil 
walk9  is  then  used  to  try  to  find  bilinear  interpolation 
coefficients  for  the  target  point  within  the  seam  sur¬ 
face  grid.  Normal  displacements  of  the  target  point 
from  the  seam  surface  grid  are  disregarded  provided 
the  target  point  is  within  a  certain  tolerance  from  the 
surface.  If  the  converged  coefficients  indicate  that  the 
target  point  is  outside  the  seam  surface  grid,  the  next 
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closest  seam  surface  grid  whose  bounding  box  contains 
the  target  point  is  tried.  This  process  terminates  when 
the  target  point  is  accepted  to  be  inside  a  seam  surface 
grid,  or  when  the  target  point  fails  to  converge  inside 
any  of  the  seam  surface  grids.  In  this  latter  case,  the 
point  is  flagged  to  be  outside  the  seam  surface  grids. 

3.3.  Creation  of  Block  Patches 

The  next  step  is  to  gather  groups  of  block  cells 
to  form  block  patches.  Let  J  and  K  be  the  index  direc¬ 
tions  for  the  marker  cells  in  a  reference  patch.  Block 
cells  from  a  reference  patch  are  related  to  each  other 
in  the  same  index  space.  On  each  reference  patch,  a 
coalescence  scheme  is  used  to  group  together  neigh¬ 
boring  block  cells  in  index  space  to  form  one  or  more 
block  patches.  Since  block  cells  from  different  refer¬ 
ence  patches  are  from  disparate  index  spaces,  a  dif¬ 
ferent  method  is  used  to  coalesce  these  cells  (see  sec¬ 
tion  3.4). 

For  each  reference  patch,  block  cells  are  first  co¬ 
alesced  in  the  J  direction  to  form  J-pencils  for  each 
K  =  constant  row  of  marker  cells  (see  Figure  4).  For 
a  particular  K  =  constant  row  of  marker  cells,  two 
block  cells  are  coalesced  if  they  are  separated  by  no 
more  than  Jc  non-block  cells.  For  example,  if  Jc  is 
zero,  then  two  block  cells  must  be  adjacent  to  each 
other  for  coalescence  to  occur.  A  larger  value  of  Jc 
would  result  in  fewer  J-pencils  at  the  end.  Figure  4 
shows  three  different  arrangements  of  block  cells  that 
would  result  in  the  same  J-pencil  after  the  J  sweep 
with  Jc  =  2.  Typically,  Jc  is  defaulted  to  a  value  of  5 
in  SBLOCK. 

After  the  coalescence  sweep  in  the  J  direction, 
the  result  is  a  collection  of  J-pencils.  Each  K  =  con¬ 
stant  row  may  contain  zero,  one  or  more  J-pencils. 
Next,  adjacent  J-pencils  in  the  K  direction  are  coa¬ 
lesced  to  form  a  block  patch  if  their  end  points  differ 
by  no  more  than  Kc  cells.  For  Kc  =  0,  the  end  points 
of  adjacent  J-pencils  have  to  match  exactly  for  coales¬ 
cence  to  occur.  A  larger  value  of  Kc  would  result  in 
fewer  block  patches  at  the  end.  Figure  4  shows  three 
different  arrangements  of  J-pencils  that  would  result  in 
the  same  block  patch  after  the  K  sweep  with  I\c  =  2. 

The  value  of  I<c  in  SBLOCK  is  defaulted  to  10  but  may 
be  modified  by  the  user.  The  boundaries  of  the  block 
patches  of  coalesced  block  cells  for  the  V-22  fuselage 
example  are  shown  in  Figure  5. 
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J  SWEEP  with  Jc=2 


K  SWEEP  with  Kc=2 


Figure  4.  Coalescence  sweep  in  J  to  form  J-pencils 
and  in  K  to  form  block  patches.  The  Jc  and  Kc  pa¬ 
rameters  can  be  used  to  reduce  the  number  of  J-pencils 
and  block  patches  created,  respectively. 


Figure  5.  Block  patches  for  the  V-22  fuselage.  Bound¬ 
aries  of  block  patches  are  marked  by  thick  black  lines. 

3.4.  Stitching  of  Block  Patches 

The  scheme  described  in  the  last  section  serves 
to  coalesce  block  cells  that  are  derived  from  the  same 
reference  patch  by  performing  connections  in  index 
space.  In  order  to  combine  block  cells  from  differ¬ 
ent  reference  patches,  connections  must  be  considered 
in  physical  space.  A  group  of  block  patches  originat¬ 
ing  from  the  same  or  different  reference  patches  are 
‘stitched5  together  to  form  a  stitched  block  patch  by 
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the  algorithms  described  below.  First,  the  candidate 
block  patches  for  stitching  are  identified.  Then,  these 
block  patches  are  joined  by  Type-1  stitching  followed 
by  Type-2  stitching.  Type-1  stitching  is  used  to  join 
pairs  of  adjacent  block  patches  that  are  almost  the 
same  size  along  the  stitch  direction.  Type-2  stitching 
is  used  to  combine  a  small  block  patch  with  an  adja¬ 
cent  large  block  patch.  These  steps  are  described  in 
more  detail  in  the  subsections  below. 


3.4.1.  Type-1  Stitching 

Type-1  stitching  is  performed  between  any  pair 
of  adjacent  block  patches  with  almost  matching  bound¬ 
aries.  For  each  point  on  each  non-periodic  boundary 
of  each  block  patch,  the  neighboring  block  patch  to 
the  point  is  identified  and  stored.  Identification  of  the 
neighboring  patch  is  accomplished  by  the  same  search 
procedure  as  the  projection  scheme  described  in  sec¬ 
tion  3.2.  Point  to  point  match  between  adjacent  block 
patches  is  not  required.  Moreover,  a  user-specified  tol¬ 
erance  parameter  is  used  to  allow  for  small  gaps  be¬ 
tween  adjacent  block  patches. 

For  the  ith  block  patch  Pi ,  let  Bi  be  a  boundary 
on  Pj,  Mi  be  the  number  of  points  on  P*,  and  My  be 
the  number  of  points  on  Bi  that  has  block  patch  Pj 
as  the  neighbor.  Type-1  stitching  will  be  performed 
between  boundary  P*  of  block  patch  P,*  and  boundary 
Bj  of  block  patch  Pj  if 

(1)  Mij/Mi  >  /,  and 

(2)  Mjj/Mj  >  f, 

where  /  is  a  threshold  parameter  set  to  0.95.  By  set¬ 
ting  /  to  be  less  than  1,  block  patches  that  almost  but 
not  exactly  match  may  be  stitched. 

The  stitching  scheme  is  designed  such  that  the 
new  patch  covers  at  least  the  same  surface  regions  as 
occupied  by  the  original  block  patches.  Moreover,  the 
grid  resolution  on  the  new  patch  is  at  least  as  fine 
as  either  of  the  original  block  patches.  The  Type-1 
stitching  steps  are  outlined  briefly  below. 

Suppose  a  Type-1  stitch  is  to  be  formed  between 
boundary  Bi  of  block  patch  Pi  and  boundary  Bj  of 
block  patch  Pj.  Let  the  number  of  points  on  P,*  and 
Bj  be  Mi  and  Mj ,  respectively,  with  Mi  >  Mj .  Points 
on  patch  Pj  in  the  direction  along  the  stitch  are  redis¬ 
tributed  uniformly  to  contain  M*  points  and  the  two 


Figure  6.  Multiple  generations  of  block  patches  cre¬ 
ated  by  Type-1  stitching.  Generation  0  (GO)  contains 
9  patches  while  generation  4  (G4)  contains  2  patches. 
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patches  are  then  concatenated.  No  redistribution  is 
necessary  if  M,  =  Mj.  Since  the  grid  points  on  the 
block  patches  are  uniformly  spaced  to  start  with,  points 
on  the  patch  boundaries  Bi  and  Bj  should  be  almost 
point  to  point  matched  after  the  point  redistribution. 
A  few  steps  of  elliptic  smoothing  is  applied  in  the  vicin¬ 
ity  of  the  stitch  to  produce  a  smooth  stitched  block 
patch.  Finally,  the  resulting  points  are  projected  back 
onto  the  reference  surface. 

After  one  round  of  Type- 1  stitching,  some  of  the 
block  patches  from  the  new  generation  may  be  com¬ 
bined  further.  Succesive  rounds  of  Type-1  stitching 
are  performed  until  no  candidates  for  Type-1  stitching 
are  found.  Figure  6  shows  9  block  patches  in  gener¬ 
ation  0.  After  4  rounds  of  Type-1  stitching,  only  2 
patches  remain. 

3.4.2.  Type-2  Stitching 

Type-2  stitching  is  performed  after  Type-1  stitch¬ 
ing  to  combine  a  small  block  patch  with  an  adjacent 
large  block  patch.  Let  Jmax*  and  ATmax*  be  the  number 
of  points  in  the  J  and  K  directions,  respectively,  for 
block  patch  i.  A  block  patch  with  Jmaxi  or  ATmaxt  less 
than  or  equal  to  7Vsman  is  flagged  to  be  ‘small’,  where 
Nsmaii  is  a  user-specified  parameter  typically  set  to  be 
about  4. 

If  J maXj  ^  Asm  all  j  then  the  J  —  1  and  J  — 
Jmaxi  boundaries  will  be  eligible  for  stitching,  while 
if  -Amaxi  ^  iV small )  then  the  K  ~  1  and  K  —  ATmaxi 
boundaries  will  be  eligible  for  stitching.  Type-2  stitch¬ 
ing  will  be  performed  between  boundary  Bi  of  a  small 
block  patch  Pt  and  boundary  Bj  of  neighboring  block 
patch  Pj  if  Mij  /Mi  >  f  where  /  is  a  threshold  param¬ 
eter  set  to  0.95  (see  Figure  7a).  The  Type-2  stitching 
steps  are  outlined  briefly  below. 

The  outer  boundary  of  the  new  patch  is  con¬ 
structed  from  outer  boundaries  of  the  original  patches 
as  shown  in  Figures  7a  and  7b.  Boundary  Bj  from  the 
large  patch  is  modified  to  enclose  boundary  B[  from 
the  small  patch.  Interior  points  from  the  large  patch 
are  retained  while  interior  points  from  the  small  patch 
are  removed.  New  points  are  inserted  in  the  new  patch 
by  linear  interpolation  in  the  interface  region  between 
the  two  old  patches  (Figure  7c).  A  few  splining  and 
smoothing  steps  are  performed  to  produce  a  smooth 
stitched  block  patch  (Figure  7d).  Points  on  the  new 
patch  are  projected  back  onto  the  reference  surface. 


Figure  7.  Steps  for  Type-2  stitching,  (a)  Original 
block  patches,  (b)  Boundaries  from  the  large  block 
patch  are  modified  to  enclose  the  small  block  patch, 
(c)  New  points  are  inserted  by  interpolation  near  the 
modified  boundaries,  (d)  Splining  and  smoothing  are 
performed  to  produce  a  new  stitched  block  patch. 
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3.5.  Grid  Spacing  Adjustments 

After  the  stitching  step,  the  surface  domain  is 
now  covered  by  a  collection  of  seam  surface  grids  and 
stitched  block  patches.  One  can  think  of  the  set  of 
stitched  block  patches  as  a  surface  domain  decom¬ 
position  of  the  regions  not  covered  by  the  seam  sur¬ 
face  grids.  The  next  step  is  to  distribute  grid  points 
appropriately  on  the  stitched  block  patches  to  form 
block  surface  grids.  Since  the  regions  of  the  surface 
domain  with  surface  discontinuities  or  high  curvature 
have  been  covered  by  the  seam  grids,  the  surface  re¬ 
gions  covered  by  the  block  grids  should  typically  be 
smooth.  Hence  a  simple  approach  would  be  to  use  a 
uniform  grid  spacing  for  the  block  surface  grids. 

Recall  that  the  marker  cells  have  been  created 
with  a  uniform  spacing  no  bigger  than  a  user-specified 
grid  spacing  A Sg.  This  spacing  has  been  carried  throu¬ 
gh  the  block  patches  and  stitched  block  patches.  One 
can  certainly  use  the  grid  points  on  the  stitched  block 
patches  as  they  are.  Then  the  block  surface  grids  are 
just  the  final  stitched  block  patches.  A  useful  enhance¬ 
ment  has  been  made  to  this  simple  scheme  by  allowing 
a  variable  grid  spacing  on  the  block  surface  grids.  Grid 
points  on  the  block  surface  grids  are  redistributed  us¬ 
ing  the  hyperbolic  tangent  function  such  that  (1)  grid 
spacings  at  the  grid  boundaries  match  the  grid  spac- 
ings  of  any  neighboring  seam  grids,  and  (2)  the  global 
grid  spacing  upper  limit  of  A Sg  is  maintained  through¬ 
out  the.  grid.  This  ensures  a  comparable  grid  spacing 
between  the  block  surface  grid  and  the  neighboring 
seam  surface  grids,  thus  allowing  good  inter-grid  com¬ 
munication  of  flow  solution  information.  Once  again, 
the  redistributed  grid  points  are  projected  back  onto 
the  reference  surface. 

3.6.  Extrapolation  to  Provide  Overlap 

The  final  step  is  to  ensure  that  the  collection  of 
seam  and  block  surface  grids  overlap  each  other.  The 
block  cells  are  constructed  such  that  they  overlap  the 
seam  surface  grids.  However,  there  is  no  overlap  be¬ 
tween  neighboring  block  surface  grids  since  they  are 
derived  from  block  cells  that  do  not  overlap  (assuming 
the  original  reference  surface  patches  do  not  overlap). 
Moreover,  small  gaps  may  sometimes  appear  between 
adjacent  block  surface  grids  from  smoothing  steps  per¬ 
formed  during  the  stitching  procedures  described  in 
section  3.4.  Excessive  smoothing  may  deform  the  grid 


boundaries  such  that  small  gaps  or  overlaps  are  intro¬ 
duced. 

In  any  case,  an  extra  layer  of  points  is  added  by 
linear  extrapolation  at  each  non-periodic  boundary  of 
each  block  surface  grid  to  provide  overlap  with  neigh¬ 
boring  surface  grids.  However,  an  extra  layer  of  points 
is  not  added  if  the  block  grid  boundary  lies  on  the  do¬ 
main  boundary  of  the  surface  geometry,  i.e.  there  is 
no  seam  surface  grid  or  other  block  surface  grid  next 
to  the  block  grid  boundary.  Finally  the  extrapolated 
points  are  projected  on  to  the  neighboring  seam  and 
block  surface  grids.  The  overlapping  seam  and  block 
surface  grids  for  the  V-22  fuselage  example  are  shown 
in  Figure  8. 


Figure  8.  Overlapping  seam  and  block  surface  grids 
for  the  V-22  tiltrotor  fuselage  (seam  grids  are  drawn 
with  thicker  and  darker  lines). 

4.  Examples 

Two  other  configurations  of  seam  and  block  sur¬ 
face  grids  are  presented  in  this  section:  a  fin-store  com¬ 
bination  and  the  X-CRV  Crew  Return  Vehicle.  The 
fin-store  case  consists  of  3  reference  surface  patches. 
The  seam  surface  grids  consist  of  a  fin-store  collar  grid 
and  a  fin  tip  cap  grid  (Figure  9a).  Six  block  surface 
grids  were  automatically  created  by  SBLOCK  (Fig¬ 
ure  9b). 

The  geometry  for  the  X-CRV  was  provided  in 
IGES  format  with  38  NURBS  surfaces  (Figure  10a). 
Although  small  gaps  exist  at  several  locations  between 
the  NURBS  patches,  no  attempt  was  made  to  fill  the 
gaps  before  proceeding  with  the  steps  below.  In  order 


777 


Copyright  ©  1997,  American  Institute  of  Aeronautics  and  Astronautics,  Inc. 


to  use  SBLOCK,  the  NURBS  surfaces  were  converted 
to  panel  networks  (PLOT3D19  multiple  surface  grid 
format)  using  NGP12  (the  National  Grid  Project  soft¬ 
ware  package).  This  step  required  only  a  small  amount 
of  effort  from  the  user.  The  next  step  involved  the  cre¬ 
ation  of  the  seam  curves  and  the  distribution  of  grid 
points  on  the  seam  curves.  There  is  currently  no  soft¬ 
ware  available  that  is  convenient  to  use  for  accomplish¬ 
ing  this  step.  Some  work  was  required  to  perform  this 
step  using  NGP. 


Figure  9.  Fin-store  configuration,  (a)  Seam  surface 
grids,  (b)  Block  surface  grids. 

Seam  surface  grids  were  then  generated  from  the 
seam  curves  using  SURGRD,  an  overset  surface  grid 
generator  on  multiple  panel  networks.9  The  7  seam 
surface  grids  required  for  the  X-CRV  were  produced 
with  very  little  effort  (Figure  10b).  It  took  about 
2  seconds  of  CPU  time  on  a  Silicon  Graphics  (SGI) 
R10000-175  MHz  workstation  to  generate  all  7  grids 
(total  of  about  8000  grid  points).  Block  surface  grids 
were  then  generated  using  SBLOCK.  A  global  maxi¬ 
mum  grid  spacing  A Sg  of  1/125  the  length  of  the  ve¬ 
hicle  was  used.  The  code  produced  32  block  patches 


after  the  coalescing  step,  and  19  overlapping  block  sur¬ 
face  grids  after  the  stitching  steps  (Figure  10c).  The 
final  19  block  surface  grids  contain  about  23000  grid 
points  and  were  created  in  about  28  CPU  seconds  on 
the  above  SGI  workstation. 


5.  Concluding  Remarks 

This  paper  describes  part  of  an  on-going  effort  to 
fill  the  void  in  software  available  to  overset  grid  users 
for  surface  domain  decomposition  and  grid  generation. 
An  algorithm  and  a  software  program  are  described  for 
performing  partial  surface  domain  decomposition  and 
surface  grid  generation  for  overset  grids.  A  complex 
surface  domain  is  covered  by  overlapping  seam  surface 
grids  and  block  surface  grids.  Given  the  surface  def¬ 
inition  and  the  seam  surface  grids,  the  block  surface 
grids  are  automatically  generated.  As  a  result,  the 
effort  required  from  the  user  is  significantly  reduced. 

The  automated  process  involves  fracturing  the 
complex  surface  domain  into  small  pieces  (marker  cells), 
and  then  using  heuristic  rules  to  coalesce/stitch  the 
pieces  back  together  into  domains  suitable  for  gener¬ 
ation  of  overset  surface  grids.  Grid  points  are  then 
automatically  placed  in  these  domains  using  simple  al¬ 
gebraic  methods.  Future  work  beyond  the  scope  of  the 
current  paper  will  include  exploring  the  development 
of  algorithms  and  graphical  user  interfaces  to  reduce 
the  effort  required  for  seam  surface  grid  generation. 


Acknowledgements 

The  authors  would  like  to  thank  Bell  Helicopter 
Textron,  Inc.  and  Boeing  Helicopters  for  furnishing 
the  V-22  geometry;  and  Ray  Gomez  of  NASA  John¬ 
son  Space  Center  for  providing  the  X-CRV  geometry. 
This  work  was  performed  as  part  of  the  DOD  High  Per¬ 
formance  Computing  Modernization  Program  funded 
through  NASA  contract  number  NAS2-14109,  Task 
23. 


778 


Figure  10.  X-CRV  Crew  Return  Vehicle,  (a)  Reference  surface  consisting  of  38  NURBS  patches,  (b)  Seam 
surface  grids,  (c)  Block  surface  grids. 
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